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We study the stochastic background of gravitational waves produced from preheating in hybrid 
inflation models. We investigate different dynamical regimes of preheating in these models and we 
compute the resulting gravity wave spectra using analytical estimates and numerical simulations. 
We discuss the dependence of the gravity wave frequencies and amplitudes on the various potential 
parameters. We find that large regions of the parameter space leads to gravity waves that may be 
observable in upcoming interferometric experiments, including Advanced LIGO, but this generally 
requires very small coupling constants. 

I. INTRODUCTION 

Gravity waves (GW) from the early universe can carry information about inflation, (p)reheating after inflation, and 
even about pre-inflation. 

During inflation, tensor modes of the classical large scale cosmological perturbations are produced from the quantum 
fluctuations of gravitons Their amplitude is proportional to the energy scale of inflation, and their spectrum 
extends over a wide range of wavelengths. Gravitational waves at the cosmological, long-wavelength part of the 
spectrum lead to B-mode polarization of the CMB anisotropy fluctuations. In high energy models of inflation (like 
chaotic inflation), where the ratio of the amplitudes of the tensor to scalar modes is about r ~ 0.1, the B-mode of 
anisotropics should be detectable by forthcoming CMB polarization experiments. 

However, in anticipation of a null signal observation of GW from inflation, one might still be able to use GW 
to constrain inflationary models in ways other than just constraining the overall energy scale. There are models of 
inflation where the total number of inflationary e-folds N exceeds the minimum required to homogenize the observable 
universe only by a small margin. For such models with anisotropic pre-inflationary expansion, gravity waves from 
pre-inflation are amplified and may contribute to the large scale CMB temperature anisotropy 0. Observational 
limits on gravity waves from CMB polarization experiments also result in constraints on N . 

If inflation occurs at lower energies (as in many hybrid inflation models), the amplitude of the resulting gravity 
waves would be too weak to be observed with CMB anistropies. However, GW from inflation at the short-wavelength 
part of the spectrum might fall in the amplitude-frequency range which accessible to future gravity wave astronomy 
projects such as BBO and DECIGO. These experiments could probe r down to the level r ~ 10~ 6 . 

Gravity waves may also carry unique information about the post-inflationary dynamics, in particular the inflaton 
decay and the subsequent evolution of its decay products towards thermal equilibrium. Often this process starts with 
preheating, a violent non-perturbative restructuring of the field configurations shortly after inflation. Preheating leads 
to large, non-linear field inhomogeneities which necessarily generate a classical GW background. This mechanism of 
GW production is complementary to the production of GW from vacuum fluctuations during inflation. In our previous 
paper Q we developed a machinery for analytic and numerical calculations of GW production in preheating models. 
There are several other papers on the subject that we will discuss below. It is also worth mentioning that the techniques 
used to study GW from preheating have parallels in calculations of GW from phase transitions / hydrodynamical 
turbulence, see e.g. Q and references therein. In 0] we developed applications of our methods to the example of 
preheating after chaotic inflation. The parameters of these large-field models are fixed by CMB normalization, which 
ensures that inflation must occur at high energy scales. As a result, the gravity waves produced from preheating in 
these models have typical frequencies of order / ~ 10 7 — 10 8 Hz, which is far beyond the observable range (see Section 
IVII for observational constraints). As conjectured in preheating after low-energy inflation may generate GW with 
lower frequencies /. One of the most popular and most studied models of inflation and preheating is hybrid inflation 
Q with tachyonic preheating Q. There are many hybrid inflation models, notably in the context of supergravity 
(see [|[ for a review) and string theory (see @ for reviews). In these models, the energy scale of inflation is not 
fixed by CMB normalization, so it was optimistically thought that for some parameters the frequency of GW from 
preheating may fall into an observable range. The subject of this paper is GW production from tachyonic preheating, 
in particular, assessing how realistic is the range of parameters that may lead to an observable signal. 

Let us briefly review the literature on the subject. The production of gravitational waves from preheating was 
originally studied in [l(| for chaotic inflation, and recently by several groups in 0, [f| [ll|, [T^, [H, Q [H, Ufl- The 
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first numerical methods [ljj, HH were based on the Weinberg formalism , which, strictly speaking is applicable 
only for isolated sources in a Minkowski background. This may lead to significant differences in the resulting gravity 
wave spectra, as explained in [3j . In [12|. the gravity wave equations in an expanding universe were solved in Fourier 
space. Another method was used in jl3l. [T3 |. where the evolution equation for metric perturbations were solved in 
configuration space. In [l3[ and the earlier version of (T3 the transverse-traceless radiative part of GW was not 
properly extracted 1 but after discussion of the problem inQ this was corrected in the second version of [l4[ . Finally, 
in we developed a method based on the Green's function solution in momentum space, to calculate, numerically 
and analytically, the production of gravity waves from a stochastic medium of scalar fields in an expanding universe. 
At present, the numerical methods of different groups [HI], 14 1 and [U seem to agree well with each other, see [HI EH 



for a comparison of the results in a model of chaotic inflation. 

Most of these papers focused on models of preheating based on parametric resonance after chaotic inflation, which 
leads to gravity wave signals at frequencies which are irrelevant observationally. Models of preheating after hyrbid 
inflation are observationally much more interesting, because they involve extra parameters and they occur at lower 
energy scales, and thus they may generate GW with lower frequencies. Gravity wave production after hybrid inflation 
was first considered in |l8j j. extrapolating the results of [Toj |. assuming that preheating occurs through parametric 
resonance. However, we know since [7[ that preheating after hybrid inflation occurs in the qualitatively different 
regime of tachyonic amplification due to the dynamical symmetry breaking, when the fields roll towards the minimum 
through the region where their effective mass is negative - a process called tachyonic preheating. It completes in one 
or very few oscillations of the fields and is very different from parametric resonance 2 . 

Before further discussion on GW, we have to comment about specific challenges for numerical simulations of 
tachyonic preheating in hybrid inflation models with several parameters 3 . These parameters include the self-coupling 
A of the symmetry breaking fields and their coupling g 2 to the inflaton. The scales of the preheating dynamics can 
be significantly different depending on the combination g 2 / A (roughly speaking, the ratio of effective masses around 
the local maximum and minimum of the potential). Therefore, care should be taken to incorporate both scales in 
the simulations. An additional subtlety is related to the initial conditions for the fields around the bifurcation point, 
namely, relatively high or low inflaton velocity, which results in qualitatively different initial regimes - quantum 
diffusion or classical fast or slow rolls around the bifurcation point, see 0] for details. The case with higher initial 
velocity is easier to model numerically. 

The production of gravitational waves properly from tachyonic preheating was first investigated in (T3 , [T3 | , for a 
specific region of the parameter space (g 2 ~ A ~ 1 and significant velocity of the inflaton at the bifurcation point), 
where the dynamical scales are of the same order. The resulting grav ity wave spectra were located at frequencies too 
high to be observable, as in chaotic inflation models. Refs. pjj, [lj| also displayed a conjecture for a low frequency 
spectrum, but the dependence on the parameters was not studied, so it remains unclear which model, if any, could 
lead to an observable signal. In 5], it was conjectured that the peak frequency and amplitude of the gravity wave 
spectrum depends in a relatively simple way on the typical scale amplified during preheating, which is usually a 
known function of the parameters. In Q we verified this conjecture numerically for a model of chaotic inflation, and 
used it to estimate analytically the peak of the gravity wave spectrum produced in two different models of preheating 
after hybrid inflation. In particular, we noticed that models with low velocity of the inflaton at the critical point are 
observationally more promising, but it is also more challeging to model this case with numerical simulations, as noted 
above. 

In the present paper, we investigate in detail gravity waves produced from preheating after hybrid inflation, focusing 
in particular on their dependence on various model parameters such as g 2 /X and the field velocities around the 
bifurcation point. We study GW production in several qualitatively different settings of tachyonic preheating, some 
of which, including the most interesting cases, were completely missed in the literature. We use a combination of 



1 Incorrect extraction of TT part, in particular, may lead to the incorrect conclusion that a significant amount of gravity waves is produced 
from the stage of scalar field "turbulence" after preheating, see Q for details. 

2 Hybrid inflation was also invoked in [l2l | to motivate a m? <f> 2 + g 2 2 x 2 model of parametric resonance at low energy scales. However, 
as noted above, this model is not relevant for preheating after hybrid inflation. The gravity wave spectra shown in 12] fall into an 
observable range for an inflaton mass of order 10 GeV, instead of 10 13 GeV in chaotic inflation (and for a coupling constant of order 
1CT 30 ). One can also consider GW production from the non-perturbative decay of a scalar field cj> in this model not necessarily related 
to the inflaton, as for example in the scenario considered in |19( . The present-day peak frequency /„ may then be estimated, in the 
same way as we do below, from the typical momentum k, ~ \/grnWoa~ 3 ^ 4 amplified in this model. If <j> dominates the energy density 
before decaying and the decay is followed by the radiation dominated era, one finds /„ > Jg 10 10 Hz, independently of the mass m 
and the initial amplitude $o of the field cf> (as long as g3> >> m for preheating to occur). This bound can be relaxed if <j> does not 
dominate the energy density, but a very small coupling g 2 is still required to achieve /* < 10 3 Hz, a necessary condition for these GW 
to be observable. 

3 Contrary to the case of resonant preheating after single-parameter chaotic inflation 
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analytical estimates and improved numerical simulations to compute the resulting gravity wave spectra and to identify 
the regions of the parameter space which may be relevant for upcoming gravity wave experiments. 

The rest of the paper is organized as follows. In Section [TTJ we briefly review the hybrid inflation models that we 
will consider. In Section [HI! we identify three qualitatively different dynamical regimes of preheating in these models 
and we make analytical estimates of the resulting gravity wave spectra. In Section ITVl we describe and advance further 
our numerical method for calculating gravity wave production from preheating. We apply this method to preheating 
after hybrid inflation in Section [V] and we study the dependence of the gravity wave amplitudes and frequencies on 
the parameters of the model. We then discuss which regions of the parameter space may lead to a detectable signal. 
Finally we discuss the implications of our results and directions for future work. Some details of our numerical 
calculations are given in the appendices. 



II. PARAMETERS OF THE HYBRID INFLATION MODEL 



We will consider simple hybrid inflation models with the potential 



V= 'x(a 2 -v 2 ) 2 



(1) 



The Higgs field a may in principle be real, complex, or consisting of any number of components, but the case of a real 
field is ruled out because it would lead to the production of dangerous sub-horizon size domain walls. We have tried 
simulations with varying numbers of a components and found the results to be largely unaffected. The simulations 
shown below are for a two-component a, for which a 2 should be understood as |<r| 2 — a\ + a\ where o\ and oi are 
two real scalar fields. We take the inflaton <f> to be real for simplicity. 



For 6 > 



where 
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(2) 



is the critical point, the fields have positive mass squared and the potential has a valley at a = 0. Inflation occurs while 
4> decreases slowly in this valley due to the uplifting term Vi n f(<f>) in (TTJ). The energy density is usually dominated 
by the false vacuum contribution, V ~ Au 4 /4. Inflation ends either at the bifurcation point when <j> = <f) c or when 
the slow- roll conditions are violated, whichever occurs first. In both cases, when <p — <p c , a acquires a tachyonic 
mass and the fields roll rapidly towards the true minimum at <f> = 0, a = v. During this rolling field fluctuations 
are exponentially excited by tachyonic preheating 0] , thus leading to a rapid decay of the homogeneous field energy, 
see Fig. [TJ This exponentially rapid growth of inhomogeneities is what drives the production of gravity waves in this 
model. 



V 




t=o 



t=80 




FIG. 1: Hybrid inflation potential. The left panel shows the potential as a function of the fields <j> and a. Only one direction 
of a is shown. The red dots (on the ridge) show the critical (bifurcation) points and the green dots (in the valleys) show the 
minima. The inflaton velocity after inflation moves it along the top of the ridge. For g 2 << A with a small initial velocity the 
fields roll along the indicated ellipse, as can be seen in the field histograms in the right panel. 
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Depending on the model, the slow-roll term Vi n f{4>) may take various forms, see e.g. Usually, it does not 

significantly affect the dynamics during preheating, except by setting the velocity with which <p reaches the critical 
point 4 . We therefore neglect Vi n j in the following and consider the initial velocity 4> c as a free parameter. We express 
this velocity in dimensionless terms as in [2fj| 

T/ _ d 4>c _ g 4>c 

Vc ~ dt ~ \v 2 W 

where a dot denotes derivative with respect to the proper time, <f> = 4>/4> c and t = mt, m = y/Xv being the natural 
mass scale for the potential. The free parameters in the model are thus A, g, v, and V c (not to be confused with the 
potential V). 

In cases where inflation lasts until <fi = <p c , it should end rapidly enough to avoid the production of Hubble scale 
inflationary perturbations of a, the so-called waterfall condition @, [2l|, [22j. Often, the whole process of preheating 
occurs in less than a Hubble time. Therefore, in the following, we will neglect the expansion of the universe. (This 
is always accurate for sufficiently small v). In this case, defining new field and spacetime variables cj> n ew = 4>/ v i 
a new = <j/v, x^^new = vx^ the value of v drops out of the field equations of motion and the amplitude of the initial 
vacuum fluctuations. Thus the independence of physical quantities is known analytically and the u-dependence of 
the gravity wave spectrum displayed below is exact. Specifically, changing v does not shift the frequencies at all and 
changes the density in gravity waves today proportionally to v 2 . 

Depending on the remaining parameters A, g and V c , the way the tachyonic instability develops in the model (|TJ) 
may occur in different regimes, as we will discuss in the next Section. 

III. ANALYTICAL ESTIMATES OF GW FREQUENCIES AND AMPLITUDES FOR DIFFERENT 

REGIMES OF TACHYONIC PREHEATING 

Before turning to the numerical calculations of GW in the following Sections, we estimate analytically how the GW 
spectra from tachyonic preheating vary with the model parameters. 

Preheating in the model (fTJ) starts when tj> = <fi c . Because of its initial velocity, the inflaton rolls classically away 
from (j) c . At the same time, the Higgs field acquires a negative mass squared which increases with time starting from 
zero at the critical point. This allows for a tachyonic amplification of its initial quantum fluctuations. When these 
fluctuations become comparable to v, the field distributions may settle rapidly around the true minimum at <f> = 0, 
a = v. The preheating process in this case has been studied in detail in [2(|, see also |23j . We will briefly review this 
case, where preheating is driven by the inflaton initial velocity, in sub-section IIII Al On the other hand, if the initial 
velocity of the inflaton is sufficiently low, the classical rolling of the inflaton may be subdominant and preheating 
may start in a different way. The initial quantum fluctuations themselves may induce a negative curvature of the 
potential around the critical point. The tachyonic amplification is then triggered by the quantum fluctuations instead 
of the inflaton's classical rolling. We will discuss this case in sub-section IIII B[ and estimate the initial velocity at 
which the dynamics crosses between these two regimes. In both cases of low and high initial velocity, when the Higgs 
fluctuations become of the order of the symmetry breaking scale, the field distributions settle rapidly around the 
true minimum if g 2 J> A. On the other hand, for g 2 << A, a significant fraction of the energy density is still in the 
relatively homogeneous inflaton, which oscillates more than once around cj> — with large amplitude. As we will see 
in sub-section IIII CI this leads to interesting differences in the process of GW production. 

Of particular interest to us will be the characteristic physical momentum fc* of the scalar fields amplified in different 
regimes of tachyonic preheating. Indeed, in Q the frequency and amplitude of the produced GW were connected to 
the dynamics of the "bubbly" inhomogeneities associated with the peaks of the random gaussian field of fluctuations 
amplified by preheating. They depend in a simple way on the typical size i?* ~ of the field bubbles. We 

verified numerically in [3j, for a model of chaotic inflation, that the main contribution to gravity wave production 
during preheating comes indeed from the violent "bubbly stage" between the linear and turbulent stages, with peak 
frequency and amplitude (in terms of the present-day GW energy density) given by 

4xl0 10 ffz 

f * * ff i/4 > ( 4 ) 
h 2 Q* gw « a x 1(T 5 (R*H p ) 2 , (5) 



4 Below we will also often call it the initial velocity at the onset of preheating, or at the bifurcation point. 



5 



where H p and p p are the Hubble parameter and the total energy density at preheating when gravity waves are 
produced. Here we have included a "fudge factor" a which may vary from one model to another (a w 0.1 for the 
X(f> 4 model considered in Q). The factor 10~ 5 arises from the redshift of the GW radiation. Not coincidentally a 
formula similar to ^ arises in the theory of GW production from the first order phase transition with the nucleation 
of bubbles of the size R. 

In many models of preheating, it is possible to estimate k*, and therefore according to (j4]) and (0, the peak 
frequency and amplitude of the GW spectrum, as a function of the parameters. We do that below for the model |T]). 
In this case, when expansion of the universe is negligible, we have 

/* ~ ^ 6 x 10 w Hz , (6) 
\v 4 

" 2n '>- ~ 2 * 10 " WW, • (7) 

where we have used p p = Xv 4 /A and H 2 = 8irp p / (3Mpj). At this level we take a w 0.1 and consider simple analytical 
estimates. The precise dependence on the parameters will be studied numerically in Section fVl 

Next, we will estimate the scale k*. It turns out that this depends, in particular, on the onset of preheating. The 
dynamics of the fields around the bifurcation (critical) point is dominated either by the classical, inertial motion of 
the field <j> superposed with the quantum fluctuations of a, or by quantum fluctuaions (quantum diffusion) of both 
fields. The estimates will depend on which process is dominant. 



A. Onset of Tachyonic Preheating by the Rolling Inflaton 

Because of its initial velocity at the critical point, the inflaton field <f> rolls classically away from <j) c along the ridge 
of the potential at a — 0. As a result, the Higgs field acquires a tachyonic mass — m 2 — Xv 2 — g 2 <fi 2 = g 2 {<f> 2 — 4> 2 ) , 
which increases in time starting from rrc£ = at the critical point. This allows for an exponential amplification of the 
quantum fluctuations of a with momenta k 2 < —m 2 ,. If the inflaton has sufficient velocity (see below), this process 
dominates the beginning of preheating. 

In this case, we can estimate the typical momenta amplified initially as follows [2^|, see also [2(| HH. Close 
to the critical point, the tachyonic mass of a increases as — m 2 — g 2 (4> 2 — <j> 2 ) — 2g 2 <fr c |0 C | At after a time At 
from the moment when <j) = (f> c . The resulting exponential growth of quantum fluctuations becomes efficient when 
At ;> y / — m 2 , leading to a typical momentum k 2 <; — mr a given by 

kl^2g 2 ^ c \ =2V c m 3 (8) 

where we have used ([3]). A broader range of momenta will be amplified as — m 2 continues to increase, but the modes 
with lower momentum ([8]) will already have an exponentially higher amplitude and their subsequent growth will occur 
exponentially faster. The process should thus be dominated by the modes with typical momentum given by l[8]). 
Inserting this result into Eqs. ©, ([7]) gives 

/, - A 1/4 V c 1/3 7 x 10 10 Hz , (9) 

h 2 n; w ~ io- 6 ^ 2/3 (^) 2 ■ (io) 

We have checked numerically that for a broad range of parameters the typical momenta amplified at the beginning 
of preheating are in very good agreement with Eq. ((5|) in the case of significant intial velocity. However, for g 2 / X << 1, 
the typical momenta are significantly shifted towards the infra-red after the first tachyonic growth. We will discuss 
this case separately below. We will see that, except in this case, our numerical results for the gravity wave spectrum 
are well described by Eqs. ([9]), (|10p. They show that, for a given non-negligible initial velocity V c , the peak frequency 
depends on the energy density Xv 4 only through A 1//4 , while the peak amplitude varies as v 2 , see They also show 
that the smaller the initial velocity V c , the lower the frequency and the higher the amplitude. However, these formulas 
cannot be extrapolated to arbitrary low initial velocity of the inflaton, as we will now discuss. 



B. Quantum Diffusion Onset of Tachyonic Preheating 



We discuss below the conditions that should apply to consider the classical inflaton rolling to be subdominant, so 
to simplify the discussion let us first consider a vanishing initial velocity at the critical point. 
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In this case, reasoning as above, one would conclude that the fields stay relatively long at the critical point = <f> c , 
(7 = 0, since the curvature of the potential vanishes there. However, quantum fluctuations of the fields around the 
critical point move the fields to the region of negative curvature slope of the potential. The resulting instability may 
be described as quantum diffusion away from the critical point of some modes among the initial fluctuations, due to 
their interactions with the other modes 0- 

In order to estimate the typical momentum amplified by this process, we first determine the direction in field space 
along which the fields are most likely to roll from the rest, i.e. the directions of steepest potential in the vicinity of 
the critical point. Letting a — vsmO and <f> = <f> c — rcosO , the potential reduces to 

V(r,0) = ^v 4 - g V\v r 3 cos8 sin 2 8 + 0{r 4 ) (11) 

at small distance r (in the fields space) from the critical point. At fixed r, this potential is minimum for tan0 = ±y/2 . 
The effective potential along these directions is given by 

V(r) = -v 4 - J-gVXvr 3 + 0{r 4 ). (12) 
4 3%/3 

Tachyonic preheating for such a cubic potential has been considered in Q- Let us briefly review their results. The 



(canonically normalised) scalar field r has initial quantum fluctuations with amplitude |£ffe| ~ l/v2fc around the 
critical point. Consider long- wavelength modes with k <J ko, for some cut-off k$. Their contribution to the mean- 
square fluctuations is given by (r 2 (k )) = (2vr)" 3 J*° d 3 k\5r k \ 2 ~ k 2 /(8-K 2 ) . Thus short-wavelength fluctuations, 
with momenta k — jko for 7 somewhat greater than one, may be seen to live on top of a quasi-homogeneous long- 
wavelength field r with an average amplitude r ~ r rms (fco) ~ ko/(2y/2n) • These short-wavelength fluctuations feel 
a negative curvature induced by the long-wavelength field r: V"(r TmB (ko)) = —Ag\T\vr vms {ko)/\/i. This may lead 
to a tachyonic amplification of the short- wavelength modes with momenta k 2 <; \V"\ ~ v2A gvko/ V3tt 2 . Taking 

for definiteness 7 p> y/2 , one may argue that fluctuations with k <, g VA v / V3ir 2 may enter a self-sustained regime 
of tachyonic growth. A more careful investigation [7[ shows that modes with somewhat higher momenta, even if 
initially more suppressed, grow faster and tend therefore to dominate the instability. We will write the characteristic 
momentum amplified by this process as 

K^CgVXv (13) 

where C is a numerical constant to be determined emprirically. 

The estimates above neglect the inflaton initial velocity. These estimates are presumably accurate when the corre- 
sponding tachyonic growth occurs faster than the one due to the classical rolling of the inflaton. Roughly speaking, 
that will be true when the typical momentum amplified is larger than the typical momentum amplified by the classical 
rolling. We thus expect the estimates above to be valid when 



^tf quant 



r 3 n 3 r 3 
* Vc<^f * <p c <^g 2 Xv 2 , (14) 



where we have used (|T5]). © and ([3]). 
Inserting (13]) into Eqs. ©, © gives 



/* - CgX 1/4 6 x 10 10 #z , (15) 

Comparing to the case of significant initial velocity we see that it is much easier to lower the gravity wave 
frequencies using small coupling constants in the case of negligible velocity V c . For g 2 = 2 A , we recover the results 
/, oc A 3 / 4 and h 2 W gw oc v 2 /X of Q. 

C. Successive GW production for g 2 << A 



We found that a qualitatively different regime of GW production from tachyonic preheating in hybrid inflation 
occurs for g 2 << A. In this case, preheating starts as in the previous sub-sections, but when the Higgs fluctuations 
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become of the order of the symmetry breaking scale, 8a ~ v, the inflaton is still relatively homogeneous, 6<f> << <p c , 
and makes more than one oscillation with large amplitude around <j) = 0. 

Indeed, the effective mass of the inflaton, = g 2 (a 2 ) is much smaller in this case, see also (22|. As a result, we 
will see that the characteristic momentum which has been amplified at the end of preheating may differ significantly 
from © or (EH). 

«p> , <|o|> <<p> , <|o|> Log [<<5<p 2 >] , Log[<<5o 2 >] 




FIG. 2: The left and middle panels show the evolution with time of the inflaton's mean normalized to (j> c (blue), (4>/4>c), and 
of the mean of the Higgs modulus normalized to its vev (red), (|<x|/w) (where \a\ 2 = a\ + erf), for X/g 2 = 0.5 (left panel) and 
X/g 2 = 2000 (middle panel). The other parameters are V c = 10" 3 and A = 10" 5 in both cases. The right panel shows the 
evolution with time of the variances Log[(S<f> 2 /cj> 2 )] (blue) and Log[(<5cr 2 /v 2 )] (red) for the case X/g 2 = 2000. 



Fig.[5]shows the evolution with time of the averages (</)) and (\a\) (in the case of significant velocity at the preheating 
onset), for g 2 ~ A (left panel) and g 2 << A (middle panel), and the same values of the other parameters. The right 
panel shows the evolution with time of the variances of the fields in the case g 2 << A. One sees clearly the large 
inflaton oscillations around the true minimum, roughly from <j) c to —<\> c for the first oscillations in the case of the 
middle panel of Fig. O Note in particular that each time the inflaton approaches ±<^ c , the minimum in the cr-direction 
is at a = 0, so that a rolls back towards the origin where the symmetry is restored. 

In both cases of Fig. [21 (|<r|) is first rapidly amplified up to a value a few times smaller than its VEV by the 
tachyonic effect as <fi rolls slowly away from <f> c due to its initial velocity. The typical momentum amplified by this 
process is given by (SJ) and is independent of g 2 /X. After that stage, non-linearities become significant. Because of 
the amplification of a, <j) acquires an effective mass which makes it roll faster towards <f> = 0. This happens much 
more slowly for g 2 << A since is much smaller in that case. As <fi decreases, the tachyonic mass of a momentarily 
increases, which amplifies further a fluctuations until backreaction shuts off the tachyonic effect. For g 2 > A, the field 
distributions then rapidly settle around the true minimum, with dispersions |<5<r| ~ v and \5<f>\ ~ 4> c . On the other 
hand, for g 2 « A, we still have \S<f>\ « <j) c , see the right panel of Fig. 

We can further study this stage by noting that the trajectory of the field distributions in field space is quite 
accurately given for some time (up to mt of order 600 in the case X/g 2 = 2000 of the middle panel of Fig. ^ by the 
ellipse 

Xa 2 +g 2 (j) 2 = Xv 2 , (17) 

which is essentially the condition dV/da = 0. This is satisfied when \cf>\ < cf> c , see Fig. [T]for illustration. Along this 
trajectory, the potential reduces to 

v = Y v2<t)2 -fx <t,4 > (18) 

and the dynamics corresponds to a single-field system with effective mass 

where (</> 2 ) ~ (0) 2 at that stage. We see that the mass squared becomes negative each time that \(j)\ > <j> c /y/3. This 
leads to the tachyonic growth of modes with typical momenta 

K~gv. (20) 

Note that this is much smaller than (5|) for g 2 « X. One can clearly see these successive tachyonic growths after the 
first, much more rapid amplification, in the right panel of Fig. [5] The successive tachyonic growths and the resulting 
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bubble inhomogeneities will lead to succesive bursts of GW productions. We will confirm this effect numerically in 
Section VC. 

In addition to these tachyonic growths, there are other non-adiabatic amplifications when m c j > (4>) ~ 0, which sig- 
nificantly affect the spectrum of 4> fluctuations. This is very similar to the preheating process in models of new 
inflation (24[. Indeed, for the potential (fT8|) the non-adiabaticity condition on the frequency, ;> ui 2 , reads 

/ \3/2 . 

Ik +m%) <J 3 |0| \4>\g i /X. This is most easily satisfied around 4> ~ ±0 c /v3, where 0. At that point, 

\(f>\ may in principle depend on the initial inflaton velocity during the first oscillations, but for the last few non- 
adiabatic amplifications a good approximation is given by \4>\ ~ \/2Xv 2 /3. This leads to the growth of modes with 
k 2 < gv around <fi ~ ±</> c /\/3. The typical momentum is lower by a factor of 2 or so, fc* ~ gv/2. 
Inserting this result into Eqs. ([6]), ([7]) gives 

K ^=A 1/4 3 x l0 10 Hz , (21) 

~ 8 xio-A(_i_)\ (22) 

Note that this does not depend anymore on the initial velocity V c . Compared to fit))) , we see that usually the peak 
frequency decreases and the peak amplitude increases for g 2 < < A . 

Since Eqs. (f2"Tj) . (f2"TJ)) do not depend on the initial tachyonic amplification, we expect them to hold also in the case 
of negligible initial velocity. Note that, in this case, it is more difficult than in (fT5|) to lower /» with small coupling 
constants. Thus it is observationally more interesting to have g 2 ;> A in the case of negligible initial velocity. 



IV. CALCULATION OF GRAVITY WAVES 

In this Section we refine the basic method of GW calculations that we developed in Q . 

The energy density of gravity waves is constructed from the amplitude of the transverse-traceless part of the metric 
perturbation h^. The TT part can be extracted by applying a momentum-space projection operator Oiji m (k) (see e.g. 
[3j), but this requires calculating the metric perturbations in momentum space. Alternatively, following [14| , we can 
solve the position-space evolution equations for the whole hij and then, when calculating gravity wave spectra, apply 
the projection operator to the end-point of the evolution results. Since the equations of motion and the projection 
operator are both linear they commute. 

The equation of motion for the metric perturbations is 



h% - V 2 h - —hij = 167rGa 3 lT£ T , (23) 

where we are using conformal time dr = dt/a and metric perturbations hij — ah^. We will use the non-TT source 
term 

Uij = -^d^drf , (24) 

and compensate by applying the projection operator (in momentum space) to the resulting hij. The other terms in 
the energy-momentum tensor vanish under the TT projection. The results of this calculation read as 

Mr,k) = ^ J\r' G[k(T,T')} a(r')T^ \r' , k) , (25) 

where Q is the Green's function of the operator d 2 + k 2 — ^— . In the typical case of power-law growth of the scale 
factor, a(t) ~ t a , the Green's function is easily constructed from Bessel functions, and may have different behaviours 
for the sub- and super-horizon modes. During preheating, the equation of state w usually jumps rapidly to a value 
close to 1/3 [25[ . In a radiation-dominated universe we can put ^- = 0. In this case the general solution (|25|) acquires 
a simple form 

16?rG 



fty(r,k) = -j— / dr' sin [fc (r - r')] a(r')T? T (t' ,k) , (26) 
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which can be matched to the source-free solution to give, at late times 

hij(r,k) = Aij (k) sin [k(r - r f )] + B tJ (k) cos [k(r - T f )] for r > r f , (27) 

where 

A tJ (k) = ^ J\r'co S [k(r f -r')]a(r')T^(r',k) 

B«(k) = ^ £ dr' sin[fc(r/-r')] a(r') ifV, k) . (28) 

Once we have Fourier transformed the metric perturbations and projected out the TT part the gravity wave energy 
density can be calculated as 

P gw = J d'kh^kyh'^k) , (29) 

where summation over i and j is understood. However, this formula will include small time oscillations of the modes, 
given by the sine and cosine terms in Eq. (|27|) above. In our Fourier space calculations we eliminated these oscillations 
by averaging over a full period, thus replacing h'^ (k)h'*j (k) with (k 2 /2) | 2 + | 2 ) . We can accomplish the same 
thing in this calculation by noting that at time r = Tf we have = hij, — h'^/k, so we can calculate the energy 
density averaged over a full oscillation as 

We use the notation |Xy | 2 = J2i j X%jX*j. 
By assuming isotropy we can write 



1 1 



dke(k*\h l3 \ 2 + \K 3 \ 2 ) ■ (3i) 



\QGa A V 

The gravity wave spectrum can then be calculated as 

4EA lf = 4^ 2| ^ |2+l ^ 12 ) ' (32) 

We make a couple of remarks about the formulas (|30)) - ([3"2"|) . As we demonstrated in 0, in the limit of t/ — > oo and 
a(t) — > const, we recover Weinberg's formula for the emission of GW from isolated sources in flat space-time. Next, 
the forms (|30|) - (|32p admit a natural interpretation in terms of a number density of emitted gravitons. Indeed, from 
the oscillating amplitudes hij (fc, r) one can construct the adiabatic invariant 

n k =^{k 2 \hj\ 2 + K\ 2 ) , (33) 

which corresponds to the number of gravitons per mode k. Then, Eq. (|30p can be re-interpreted in terms of the energy 
density carried out by the gravitons 

P - = 32^4 i dZkkn *- ( 34 ) 



k n k 



Correspondingly, Sk — s( . ; , 

From Eq. (|32p . the derivation of today's spectrum occurs as it did with our Fourier space calculation, adapted for 



the hybrid inflation case as discussed below. In the chaotic inflation case we evolved the spectrum to today by defining 
a time tj at the end of the simulation after which we took the equation of state to be w — 1/3. In the case of hybrid 
inflation, we are not considering expansion during preheating and we assume that the equation of state approaches 
w = 1/3 within a few Hubble times afterwards, so we effectively take tj to be the end of inflation, and simply evolve 
to today's variables with 

f=A I (4xlO w Hz), (35) 
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fl qw h 2 = (9.3 x 10- 6 )^^ , (36) 
9 

where we take p = (1/4) Aw 4 . 

In the next Section we will describe results of numerical calculations of GW from lattice simulations of tachyonic 
preheating. However, the infra-red (IR) tail of the spectrum, which is often interesting for observations, is not captured 
by the finite-size simulations. Eq. (|26[) is useful for deriving the IR asymptotics of the gravity wave spectrum. For 
momenta small compared to the peak of the scalar field spectra, k << fc», the convolution of the scalar field spatial 
derivatives is independent of k, T^ T (r',k) ~ T^ T (r',0). This corresponds to the quadrupole approximation. We 
can then distinguish two different regimes, depending on the time variation of the Green's function sin [k (t — r')] 
and the source Ty T (r',0). For sufficiently small k, the Green's function varies more slowly than the source and may 
be taken out of the integral in (f2"6")l . This gives fl gw oc f 3 for the IR tail of the GW spectrum (f3"6")l . (13"2")) . On the 
other hand, there can be an intermediate frequency range between the IR tail and the peak where sin [k (t — r')] 
varies more rapidly than T? t (t',0). The time integral in ([26")) then gives an extra 1/k factor, leading to Q, gw oc / 
for some frequency range below the peak. In the transitional region we expect a power-law slope, £l gw oc / 7 with 7 
varying between 1 and 3. For the model of parametric resonance that we considered in 0], we observed the £l gw oc / 
behaviour. For tachyonic preheating, however, the source varies typically with a characteristic time given by 1/fc*, 
i.e. more rapidly than the Green's function for any k < k*. Therefore, we expect the IR tail Q gw oc / 3 to be valid 
quickly below the peak. For the IR part of the spectrum just below the peak that we can probe in our simulations, 
we found 7^2 — 2.5. The exact value was different for the three cases of significant initial velocity, negligible initial 
velocity and g 2 « A, but it was the same in all our numerical results for a given case. 

V. NUMERICAL RESULTS AND PARAMETER DEPENDENCE 

In this section we report the results of numerical calculations of GW radiation from tachyonic preheating after 
hybrid inflation. Some technical details of the numerical simulations are described in the Appendices. In particular, 
we study how the gravity wave characteristics depend on the parameters of the hybrid model. This will allow us to 
directly check our analytical estimates and find fits for the peak frequency and amplitude. The analytical estimates 
of Section IIIII were based on the picture of a sudden growth of the density "bubbles" associated with the high peaks 
of the random gaussian fields of initial quantum fluctuations. Therefore we begin this Section with an illustration of 
the growth of these bubbles, shown in Fig. (3) The four panels show the field amplitude (left of each panel) and GW 
energy density (right of each panel) at four moments of time on a two-dimensional slice through the three-dimensional 
lattice. The initial bump of the field produces an exponentially growing bubble, and we can clearly see a burst of GW 
radiating from this growing bubble. The background stochastic GW radiation is the superposition of such concentric 
bursts. 

The configuration space picture of Fig. [3] is complemented by the momentum space spectra of gravity waves. Fig. [4] 
shows the gravity wave spectra from the lattice simulation with one representative set of parameters corresponding 
to the case A/5 2 = 1/2 and the onset of preheating dominated by classical rolling. All spectra shown in this paper 
are scaled to present-day units. The different curves show cumulative results for different moments of times during 
the simulation. For some time after inflation no gravity waves are produced, then there is a burst of growth during 
preheating, and then the spectrum saturates at a stationary level. A clear peak frequency is established early in the 
growth and remains roughly constant as the amplitude grows. 

As noted above, the model ([T]) involves essentially four parameters related to preheating: the coupling constants A 
and g, the symmetry breaking VEV v, and the unitless field velocity at the critical point, V c . 

First, we discuss how GW spectra calculated from lattice simulations depend on the parameters in the case where 
A = 2g 2 . From the point of view of numerical simulations this is the simplest case because there is no gap between 
different mass scales. Later we will study how the results depend on the ratio g 2 /X- 

As noted in Section [TH the dependence of the spectra on v can be calculated analytically in the absence of expan- 
sion. Specifically, changing v does not shift the frequencies at all and changes the density in gravity waves today 
proportionally to v 2 . This result should hold as long as v is small enough to maintain the waterfall condition. 

As we anticipated in Section lTiTl variation with V c , the inflaton velocity at the bifurcation point, results in significant 
variations of the GW spactra. For the case with significant initial velocity Eqs. ©, (|10j) accurately predict the peak 
frequency and amplitude to within an order of magnitude for all of the parameters we have tested, which includes 
initial velocities ranging from V c — 10 -5 to 10~ 2 for A = 10~ 14 and 10 -5 . 

Fig. [5] shows the final gravity wave spectra for a range of initial velocities for A = 10~ 14 . Increases in V c lead to 

increases in peak frequency and decreases in peak amplitude, but the shape of the spectra are quite similar. Fig. [5] 

1/3 

shows the numerical peak frequencies for the cases from Fig. [S] Eq. © predicts that f pea k oc V c so the figure shows 
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FIG. 3: Field amplitude and gravity wave density as a function of space on a two dimensional slice through the lattice. The slice 
is chosen at the height where the first Higgs bubble appeared. The simulation is of the model ([TJ with A = 10 -5 , X/g 2 = 0.5, 
v = 10" 3 , and V c = 0. 

Log(n gw /r) 





Log© (Hz) 

6.0 6.5 7.0 



FIG. 4: Gravity wave spectra for the model JT} with A = 10~ 14 , X/g 2 = 0.5, v = 10 -3 , and V c = 10" 5 . Lower (red) curves 
correspond to earlier moments of the simulation and higher (bluer) curves are from later moments of the simulation. All results 
are scaled to the present-day spectrum. 

a best-fit line with slope 1/3, illustrating that the dependence matches the prediction. The actual values of the peak 
frequencies are all within a factor of 3 of the Eq. ©, which is in any case only an order of magnitude estimate. 

Eq. © is not expected to be accurate in the limit of very low initial velocities, as discussed in sub-section IIII Bl 
For example, Figs. [7][8]show the spectra and peak frequencies for A = 1CP 5 . For this larger value of A the formula 
(J9j) fails for V c <; 10~ 4 . Figure [7] shows that the spectra for V c = 10 -5 and V c — 10~ 6 are nearly identical to the 
V c = spectrum, so below V c < 10~ 4 there is no V c dependence. The fit shown in this figure once again has the slope 
predicted by Eq. ([9]) and a height adjusted to fit the data, but in this case the fit is to the last two points only. The 
actual values for peak frequency for these last two points match Eq. (J9j> to within a factor of 2. 

For lower A, such as the results shown above, the initial velocity for which quantum diffusion dominates is much 
lower and thus the results in Fig. [6] show the V c dependence discussed above. 

For g 2 = 2A, our estimate (fT4")) indicates that the velocity V c could be considered negligible below the cutoff 
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FIG. 5: Dependence of GW final spectra on V c . The 
curves from left to right correspond to V c = 10 -5 , V c = 
10" 4 , Vc = 10" 3 , and Vc = 10" 2 . Here A = 10" 14 . 
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FIG. 6: The peak frequency for the gravity wave spec- 
trum as a function of initial velocity for A = 10 _1 . The 
line shows the slope of the prediction Q with the overall 
height adjusted to fit the data. 



J 8.4 - 




FIG. 7: Final gravity wave spectra after preheating 
for A = 10" °. The plots from left to right show 
Vc =: 0, 10" 6 ,10" 5 ,V C = 10" 4 ,10" 3 ,10" 2 . The plots 
for V c = 0, 10 -6 , and 10 -5 lie almost perfectly on top of 
each other. 



FIG. 8: The peak frequency for the gravity wave spec- 
trum as a function of initial velocity for A = 10 . The 
line shows the slope of the prediction Q with the overall 
height adjusted to fit the last two points. The horizontal 
line is the peak frequency for V c = 0. For V c < 10 -4 the 
peak remains constant at the V c — value. 



V c Rj C 3 A 3 / 2 . From the data above we conclude that for A = 10 -5 this cutoff occurs between V c = 10~ 5 and 
V c = 10 , which means C ~ 10. It is hard to estimate, however, if this value is independent of A. For lower values 
of A we were unable to simulate slow enough initial velocities to see the crossover between these two regimes. We can 
thus only conjecture that this value C ~ 10 will be roughly constant for different values of A. This conjecture must 
be true if our predicted dependence f pea k oc A 3 / 4 is correct for the case of negligible initial velocity. 

Next, we investigate an impact of variation with A on the GW production. For the case with significant initial 
velocity V c the equations above predict that changing A should not change the amplitude but should simply shift the 
frequency as f pea k oc A 1 / 4 . Fig. [5] clearly shows this to be the case. 

As discussed above, we can not numerically test the dependence of the spectrum on A for the case of negligible 
initial velocity, but we believe our predictions should be accurate for this case because at low initial velocity the 
bubble description of preheating is accurate. 

We now discuss how our numerical results for the gravity wave spectrum depend on the ratio g 2 /\. For a given 
inflationary potential Vi n f(</>), changing g 2 /X changes the critical point <j) c and the unitless velocity V c at that point, 
see ([3]). However, in order to isolate the g 2 /X dependence, we will keep V c constant in this section. We first consider 
the case of non-negligible initial velocity. 

In this case, gravity wave spectra for different values of g 2 /X are shown in Fig. [TD1 for V c = 10~ 3 . For g 2 > A, the 
peak frequency and amplitude are roughly independent of g 2 /A, in agreement with our original prediction, Eqs. ([9]), 
([ID)). (When g 2 increases, rescattering effects become more important and the UV part of the GW spectrum has 
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FIG. 9: Gravity wave spectra for A = 1CP 14 (left) and A = 10~ 5 (right), each with V c = 10 -3 . The vertical line through the 
A = 10 -5 peak marks the maximum of that peak. The vertical line through the A = 10 -14 peak is a prediction based on 
extrapolating the A = 10 -5 peak frequency assuming f pea k oc A 1//4 . 

higher amplitude). On the other hand, for g 2 << A, it is clear from Fig. [TO] that decreasing g 2 /X decreases the peak 
frequency and increases the peak amplitude, in aggreement with Eqs. (|2ip . ([22]) . Note also that, for a given value of 
g 2 / X << 1, the dependence on A is very well given by /* oc A 1 / 4 and h 2 f2* w independent of A, again in agreement 
with Eqs. (|2"Tj) . (|2"2"]l . Note that Fig. [TO] shows the gravity wave spectrum for v = 10~ 3 Mpi. For g 2 « A, a lower 
value of v may be necessary to consistently neglect the expansion of the universe, as discussed below. Lowering v 
leaves the shape of the spectra unchanged but reduces their amplitude proportionally to v 2 . 
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FIG. 10: Numerical results for the gravity wave spectra for several values of X/g 2 , for A = 10~ 14 (left) and A = 1CP 5 (right), 
each with V c = 10" 3 . The spectra for A = 10" 5 (right) are, from top to bottom, for X/g 2 = 20000 (black), 5000 (blue), 500 
(red), 50 (black), 0.5 (blue), 0.005 (red) and 0.0005 (black). The spectra for A = 10" 14 (left) are for X/g 2 = 5000 (top, blue) 
and 500 (bottom, red). To ease the comparison, all the spectra are shown for v = 10 Afpi, although for g 2 « A a lower 
value of v may be necessary to consistently neglect expansion of the universe (see the main text for details). This lowers the 
spectra as h 2 il gw oc v 2 . In other words, the {/-axis of the plot is really the log of h 2 fi gw /(10 3 v/Mp\) 2 . 

Decreasing g 2 / X results in qualitatively different GW production. Fig. [TT] shows the accumulation with time of 
the total energy density in gravity waves for two cases, the left panel for g 2 ~ A and the right panel for g 2 << X. 
In the first case we see a single burst of GW production from tachyonic preheating. In the second case, we clearly 
see successive bursts of gravity wave production, due to the successive tachyonic and non-adiabatic amplifications 
discussed in subsection IIII CI and the subsequent bubble collisions. Because the characteristic time scale in this case 
is given by mvX/g, for g 2 << A preheating and gravity wave production take more time than for g 2 ~ A. As a 
result, the upper bound on v necessary for the expansion of the universe to be negligible is lower than for g 2 J> A. 
This decreases the maximum GW amplitude as h 2 tt gw oc v 2 . If the peak of the gravity wave spectrum is reached 
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after a proper time r , we require it to be much smaller than the Hubble time 
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FIG. 11: Evolution with time of the total energy density in gravity waves, p gw /(10 3 v/M-p\) 2 , for X/g 2 — 0.5 (left) and 
X/g 2 = 2000. The parameters are the same as in Fig. [2] 

In fact, for all the cases with g 2 << A that we considered, the peak amplitude of the GW spectrum continued to 
slightly but constantly increase with time after the last burst of GW production (i.e. after mt ~ 1000 for the right 
panel of Fig. Hip , and the peak frequency tended to move further towards the infrared. However, if we consider very 
late times expansion will become significant, thus invalidating the results of our simulations. For low values of v this 
growth could continue longer without being diluted by expansion, but the overall amplitude would be lower. For each 
set of parameters A, g, and V c there is thus an optimal value of v that will produce the greatest amplitude of gravity 
waves before expansion becomes significant. For each set of parameters A, g, and V c we define v op t to be this optimal 
value and we define the time r as the time at which expansion would become significant for v = v op t- The spectra 
shown in Fig. [TOl were obtained at that time. 
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FIG. 12: Peak frequency /* (left), peak amplitude h 2 fi| w (middle) and time mr defined in the main text (right), for V c = 10 3 , 
A = 10~ 5 , v = 10 -3 , and several values of g 2 /X << 1. The lines shown are the best linear fits for all the points except the one 
with X/g 2 = 20000 (extreme right): f, oc (X/g 2 )~ 1/2 , h 2 fl* g „ oc {X/g 2 ) 1 - 16 and mr oc {X/g 2 ) ' 53 . 

Fig. [T2] shows in the left panel: the peak frequency, in the middle: the peak amplitude and in the right: the 
time mr defined above, for V c = 10~ 3 , A = 10~ 5 and several values of g 2 /X << 1. We see that the prediction 
([2~Tj). /» cx g, is very well satisfied. The peak amplitude grows slightly more rapidly than (f2"2"| when g 2 /X decreases, 
h 2 cx (X/g 2 ) 1 16 . Finally, mr grows slightly more rapidly than y/X/g, mr oc (X/g 2 ) ' 53 . In these fits, we 
disregarded the point at the extreme right (X/g 2 = 20000), which is slightly displaced compared to the others. In fact, 
for X/g 2 J> 10 4 , and for the initial velocity V c — 10~ 3 considered here, the inflaton reached \<j>\ > 4> c and spent most 
of the time in these regions. Indeed, in these regions, the minimum in the cr-direction is at a — and the effective 
mass for (f> is very low, m 2 ^ — g 2 (<J 2 ). Apparently, this slightly increases the peak amplitude and decreases the peak 
frequency. However, it may not be a good approximation to neglect Vi n f(4>) in the region where \<f>\ > 4> c ■ Also, we 
decreased g 2 /X while keeping V c fixed. However, for a given physical initial velocity 4> c , decreasing g 2 /X decreases 
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the unitless velocity V c , see ([3]). For sufficiently small V c , the inflaton will not reach values \4>\ > (f> c . We will thus 
neglect this effect below, although it may be useful to further decrease the peak frequency. 

We then obtain the following fits for the peak frequency and amplitude of the GW spectrum for the case g 2 /X << 1 

/* ~ 4= A 1 / 4 10 10 - 25 Hz (38) 
v A 



A 

? J \M P1 



h 2 figw - 10- 5 - 5 -5 ) ( — ) ■ (39) 



From the fit of m r, the constraint (f3"T|) for the expansion of the universe to be negligible becomes 



A ^ - a53 



<< kt 1 - 3 . (40) 



Mpi V.9 

Note that this condition always implies (j> c << Mp\ . Together with (|3T)|), it also implies h 2 fl* w << 10~ 81 (A/5 2 ) 01 . 

For instance, for A ~ 0.1, /* <; 10 3 Hz for g 2 <^ 10~ 15 , and the amplitude ([39]) may be relevant for Advanced LIGO 
while satisfying the constraint (|40|>. 

We expect similar results as above in the case of negligible velocity V c . In this case, where quantum diffusion 
dominates, observationally more interesting results can be obtained without the requiremenet g 2 << A. 



VI. DISCUSSION AND PERSPECTIVES 



In this paper, we studied the stochastic background of gravitational waves produced from tachonic preheating after 
hybrid inflation in the very early universe. The present-day frequencies and amplitudes of these gravity waves may 
cover a wide range of values, depending on the main phenomenological parameters introduced in Section II: the VeV v 
and self-coupling A of the symmetry breaking fields, their coupling to the inflaton g 2 and the (unitless) initial velocity 
of the inflaton at the critical point V c . We developed analytical and numerical tools to calculate the resulting GW 
spectra and to study in detail how they depend on these parameters. We identified three dynamical regimes leading 
to qualitatively different results for the GW spectra produced from tachyonic preheating: (i) the case with g 2 A 
and the onset of preheating driven by the classical rolling of the inflaton, (ii) the case with g 2 > A and quantum 
diffusion at the onset of preheating (corresponding to a negligible initial velocity of the inflaton at the bifurcation 
point) and (iii) the case with g 2 << A where GW are produced in successive bursts. Only the first case was considered 
in previous works about GW from tachyonic preheating [T3 . [Tij , where the resulting GW spectrum was computed in 
the specific range g 2 ~ A ~ V c ~ 0(1). 

Based on our results, we can determine the range of parameters of hybrid inflation models for which preheating may 
lead to a GW signal that is potentially observable. Let us first discuss the frequencies and amplitudes of stochastic 
GW backgrounds that are relevant for GW astronomy. We consider the present-day frequencies / and the spectrum 
of energy density per logarithmic frequency interval h 2 f2 gw (/). Fig. 1131 shows the sketch of expected sensitivities of 
planned and future interferometric experiments: Advanced LIGO, LISA, the Einstein Telescope, Big Bang Observer 
and DECIGO (see for the relevant references). Note also the recent atomic interferometric sensor proposal of [2^ |. 
which may be sensitive to GW backgrounds with frequencies in the range 1 — 10 Hz and in the LISA range. We also 
show the BBN and ms pulsar bounds, and predictions for the stochastic GW background generated from inflation, 
for different values of the parameter r (the ratio of the tensor and scalar amplitudes of the inflationary cosmological 
perturbations). Cosmological GW signals will be obscured by several expected astrophysical foregrounds, notably 
from White Dwarf binaries [28[ , see Fig. [T3J 

There are also several bar and spherical resonant detectors (see e.g. [29| ') operating in the kHz range. Other 
experiments have been proposed at higher frequencies, up to 100 MHz [30]. However, the sensitivity to h 2 f2 gw 5 drops 
dramatically when the frequency increases. Indeed, we also have to take into account the current "Standard Quantum 
Limit" 



5 The concept of h 2 f2 gw for stochastic backgrounds is not applicable for GW signals peaked within a very narrow frequency band. 
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FIG. 13: Expected sensitivities of interferometric experiments (adapted from [26j]) compared to gravity wave spectra from 
tachyonic preheating in the cases of significant intial velocity with g 2 > A (green), negligible initial velocity with g 2 > A (red), 

and g << A (blue). Also shown are the BBN and ms pulsar bounds, the stochastic background from extragalactic White 
Dwarfs binaries (taken from HH), the inflationary background for two values of the tensor to scalar ratio, and the spectrum 
from preheating in a A (f> 4 model of chaotic inflation (double spectra in grey at the right, calculated with our previous code @|). 
The GW spectra from tachyonic preheating, for the peak from left to right, correspond to the following values of the parameters: 
(1) A = 2g 2 — 10" 14 , v = 3 10" 7 M PJ and negligible initial velocity ; (2) A = 0.1, g 2 = 1(T 16 and v = 3 10~ 10 M PJ (independent 
of V c ) ; (3) A = 2g 2 = 10" 14 , v = 10" 3 M vx and V c = 10 -3 ; (4) A = 0.1, g 2 = 10" 4 and v = 6 10" 4 M Pl (independent of V c ) ; 
(5) A = 2g 2 — 10" 5 , v = 10" 3 Mpi and negligible initial velocity. 

shown in Fig. 1131 The region on the right of this line is not observable with interferometric experiments due to the 
shot noise fluctuations of photons. Here x is the squeezing parameter, which experimentalists are trying to push from 
one to a few. 

Exemples of GW spectra from tachonic preheating calculated in this paper are shown in Fig. [T2J We also show 
in this figure a GW spectrum from preheating after chaotic inflation calculated in [3j (the double-spectrum in grey), 
which is located in the high-frequency side and is not observable. GW from tachyonic preheating may cover a much 
wider range of frequencies, spreading from the highest frequencies shown in the figure to the region in between the 
SQL limit and the WD binaries, which is in principle observable 6 . In each case shown in the figure, a lower Higgs 
VeV would lead to lower amplitudes and smaller coupling constants would lead to smaller frequencies. 

Our analytical and numerical results for the peak frequency /* and the peak amplitude h 2 f2* w of the GW spectra 
are well described by Eqs. (f5|- (fl"U|) in the case with g 2 J> A and the onset of preheating driven by the classical rolling 
of the inflaton, Eqs. (|T5|) - (p^j) in the case with g 2 ;> A and quantum diffusion at the onset of preheating (negligible 
initial velocity of the inflaton at the critical point) and Eqs. (f3"5|) -([3T) l) in the case with g 2 << A. The peak frequency 
depends essentially on the coupling constants and is independent of the Higgs VeV v. As a result, A ~ g 2 ~ 0(1) 



For the sake of curiosity, one can notice the analogy with the 2.73K CMB radiation which is observable in the strip of frequencies 
between the borders of the foreground dust and synchrotron radiations, as was outlined in early 1960. 
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leads to GW at very high frequencies, independent of the energy scale during inflation, and the only way to lower 
the frequency is to lower the coupling constants. On the other hand, h? f2* w oc v 2 , so that at any frequency the 
amplitude can be relatively high for sufficiently high v, roughly up to the maximal bound 7 h 2 fl gw < 10~ 6 . Such a 
high energy density in GW implies in particular that these GW may already be observable by Advanced LIGO, but 
this generally requires very small coupling constant(s), see Fig. Q2J For illustration, we show in Fig. Q3J for each of 
the three cases discussed above, the range of parameters of the hybrid inflation models such that the peak frequency 
of the GW produced from preheating satisfy /* < 10 3 Hz. This is basically the condition for these gravity waves to 
be observable if v is set to the maximum possible value consistent with the waterfall condition. In the first case of 
Fig. 1141 corresponding to g 2 ^> A and significant initial velocity, the coupling of interest A typically has to be extremely 
small, down to A < 10 -30 for V c ~ 1. In the second case, with negligible initial velocity, we may have g 2 ~ A ~ 10 -11 . 
In the third case (g 2 « A), the most interesting regime corresponds to A ~ 1, but this still requires g 2 < 10~ 14 . 




-10 -8 -6 -4 -2 -10 -8 -6 -4 -2 2 4 6 8 10 12 14 

U>g[VJ Log[A/ s 2 | Logli/S 2 ] 



FIG. 14: The regions of parameter space for which the peak of the gravity wave signal satisfies /» < 10 3 Hz. The panels show: 
Left - significant initial velocity and g 2 > A. Middle - negligible initial velocity and g 2 > A. Right - j 2 < A. 

Above we focused on the peak frequency /* of the GW spectra. At lower frequencies, the spectra have a power-law 
tail, £! gw oc f a . For tachyonic preheating, we found in this paper an intermediate regime with a relatively steep slope, 
a = 2 — 2.5, before the IR tail with a = 3 corresponding to the modes produced outside the Hubble radius. The 
infra-red part of the spectrum does not significantly improve the constraints on the parameters shown in Fig. 1141 As 
outlined above, the range of parameters leading to the maximum amplitude in GW correspond to fc* slightly below 
H. In this case, the spectra decrease as fi gw oc / 3 for frequencies just below the peak. In particular, since this is the 
same dependence as the SQL limit ([4"Tj) . if the peak of the GW spectrum is located at the right of this line, the same 
is true for the IR tail. 

In sum, our sober conclusion is that GW from tachyonic preheating after hybrid inflation can fall in the observable 

range of the h 2 il gw / plane that is currently conceivable only for uncomfortably small values of the coupling 

constants of the model. 

We should note, however, that particle physics models of hybrid inflation with such small coupling constants have 
been considered. First, if the field a in the model |T]) corresponds to a flat direction, we may expect its negative mass 
squared at the origin to be of the order Aw 2 ~ TeV 2 , and the true minimum of the potential to be located at a VeV 
at a much higher scale, as high as v ~ Mpj. This gives a coupling constant as low as A ~ 1CP 30 . Hybrid inflation 
models along this line have been proposed in (3ll. J32I. Another model which, for the same reason, naturally involves 
very small coupling constants is thermal inflation [331 ] . Preheating mechanisms in this model have not been studied in 
detail yet, but tachyonic amplifications may lead to an interesting GW signal, see also [34|]. Another model of hybrid 
inflation in supergravity is P-term inflation (35| , which includes F-term and D-term models as special cases and may 
also be realized in D3/D7 models in string theory. Satisfying the observational bounds on cosmic strings in P-term 



7 This bound follows from the approximate equation (|5j with the constraint fc* > H (meaning that the scalar fields and the gravity waves 
are produced inside the Hubble radius). The same constraint, fe, > H, is also roughly the condition for the expansion of the universe 
to be negligible and for the waterfall condition to be satisfied. It also implies that v « Mpj and (for the small coupling constant we 
consider) that the inflationary model is not ruled out by the non-observation of inflationary gravity waves. 
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inflation may require small couplings (A ~ g 2 < 10~ 10 ), see [33. Finally, we note that the case g 2 << A may have 
connections with brane-anti-brane inflation in a warped throat |36| , which is another prototype of hybrid inflation in 
string theory. According to [37j j. the end of inflation in this scenario should exhibit similarities with the model ([1]) 
with g 2 << A (see in particular Eqs. (30)-(32) of the second paper in [ST]])- In this case, a very small g 2 would emerge 
naturally because of the exponential warping of the throat. 

More generally, the hybrid inflation model {!]) that we considered in this paper is only one of the much broader 
class of inflationary models where preheating occurs through tachyonic effects. GW production from other models 
involving tachyonic preheating is a subject for further investigations. 
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Appendix A: Numerical Calculations 



Because of the range of scales needed to accurately determine the spectra for our runs the simulations shown here 
were done using CLUSTEREASY |H, the parallel programming version of LATTICEEASY ^|. We evolve the 
non-TT metric perturbations along with the scalar fields using Eq. (|23[) and then extract the TT part and calculate 
the gravity wave spectrum as describe in Section IIVI 

For runs with significant initial velocity, we did some runs with the initial conditions described in [2(| , but we found 
the results virtually identical to those obtained by simply starting with vacuum initial conditions at the critical point, 
so the simulations shown here were all done starting from vacuum modes. 

We found that the gravity wave spectrum is highly sensitive to the UV and IR tails. When the lattice has insufficient 
UV coverage this results in a large, spurious growth in the UV part of the spectrum. This growth can in turn raise 
or lower the rest of the spectrum, leading to a generally unreliable spectrum. Likewise, having insufficient IR, even 
well below the region of the peak, can lead to inaccurate results in the region of the peak. These effects are for the 
most part not a result of the initial spectrum. For most of the physical parameters described here we did runs with 
a wide range of cutoffs in the initial spectrum and found the final results virtually independent of these cutoffs. In 
short an accurate spectrum simply required a grid with enough modes in a wide band around the peak. 



Appendix B: Symmetries of the Metric Perturbations 



In principle the energy density in gravity waves (|30[) involves a sum over the nine components of , but this can be 
reduced through symmetry. The tensor is symmetric (hij — hji, 3 equations), traceless (ha — 0, 1 equation), and 
transverse (fcj/iy = 0, 3 equations), which reduces the number of independent components from nine to two. However, 
the two components that you need to find the rest are different for different k because some of components are zero 
along certain directions in k space. The scheme we used for calculating fl gw for all directions in k is the following: 

For all points with k z ^ (off the k x , k y plane) calculate h\\ and /112 and the symmetry equations give 

hijh* 3 = 2 1 2 3 [(fc 2 + kl) h\ 1R + 2kik 2 hii iR h 12 ,R + (fcf + fc 2 ) h\ 2R ] + equivalent terms for h ijt i (42) 

For all points with k z = and k y ^ (on the k x , k y plane but off the k x axis) calculate hu and and the 
corresponding solution is 

k 2 + k 2 

2 , 4 [(k 2 + fcf) h 2 x R + k 2 h 2 3 r] + equivalent terms for h^j (43) 
k 2 

For all points with k y — k z — (the k x axis) calculate /122 and /123 and the corresponding solution is 

2 (^22 r + ^-23 r) + equivalent terms for j (44) 
Analogous equations are used to calculation h'^h'^. 



19 



[1] A. A. Starobinsky, JETP Lett. 30, 682 (1979) [P isma Zh. Eksp. T eor. Fiz. 30, 719 (1979)]. 
[2] A. E. Gumrukcuoglu, L. Kofman and M. Peloso. larXiv:0807.133"5l [astro- ph], 

[3] J. F. Dufaux, A. Bergman, G. N. Felder, L. Kofman and J. P. Uzan, "Theory and Numerics of Gravitational Waves from 
Preheating after Inflation," Phys. Rev. D 76, 123517 (2007) arXiv:0707.0875l [astro-ph]]. 

[4] T. Kahniashvili, A. Kosowsky, G. Gogoberidze and Y. Maravin, Phys. Rev. D 78, 043003 (2008) (arXiv:0806.0293l [astro- 
ph]]. 

[5] G. N. Felder and L. Kofman, "Nonlinear inflaton fragmentation after preheating," Phys. Rev. D 75, 043518 (2007) 
[arXiv:hep-ph/0606256j . 

[6] A. D. Linde, "Axions in inflationary cosmology," Phys. Lett. B 259, 38 (1991) ; A. D. Linde, "Hybrid inflation," Phys. 

Rev. D 49, 748 (1994) [arXiv:astro-ph /9307002] . 
[7] G. N. Felder, J. Garcia-Bellido, P. B. Greene, L. Kofman, A. D. Linde and I. Tkachev, "Dynamics of symmetry breaking 

and tachyonic preheating," Phys. Rev. Lett. 87, 011601 (2001) [arXiv:hep-ph/0012142j ; G. N. Felder, L. Kofman and 

A. D. Linde, "T achyonic instability and dynamics of spontaneous symmetry breaking," Phys. Rev. D 64, 123517 (2001) 

[arXiv:hep-th/0106179| . 

[8] D. H. Lyth and A. Riotto, "Particle physics models of inflation and the cosmological density perturbation," Phys. Rept. 

314, 1 (1999) [arXiv:hep-ph/9807278j . 
[9] F. Quevedo, Class. Qua nt. Grav. 19, 5721 (2002) [arXiv:hep-th/0210292| . R. Kallosh, Lect. Notes Phys. 738, 119 (2008) 

larXiv:hep-th/0702059| . C. P. Burgess, PoS P2GC, 008 (2006) [Class. Q uant. Grav. 24, S 795 (2007)] |arXiv:0708. 2865 

[hep-thjj. L McAllister and E. Silverstein, Gen. Rel. Grav. 40, 565 (2008) ]arXiv:0710.2l)5Tl [hep-th]] . 
[10] S. Y. Khlebnikov and I. I. Tkachev, "Relic gravitational waves produced after preheating," Phys. Rev. D 56, 653 (1997) 

[arXiv:hep-ph/9701423] . 

[11] R. Easther and E. A. Lim, "Stochastic gravitational wave production after inflation," JCAP 0604, 010 (2006) 
[arXiv:astro-ph/0601617| . 

[12] R. Easther, J. T. . Giblin and E. A. Lim , "Gravitational Wave Production At The End Of Inflation," Phys. Rev. Lett. 99, 

221301 (2007) [arXiv:astro-ph/0612294] . 
[13] J. Garcia-Bellido and D. G. Figueroa, "A stochastic background of gravitational waves from hybrid preheating," Phys. 

Rev. Lett. 98, 061302 (2007) [arXiv:astro-ph/0701014] . 
[14] J. Garcia-Bellido, D. G. Figueroa and A. Sastre, "A Gravitational Wave Background from Reheating after Hybrid Inflation," 

Phys. Rev. D 77, 043517 (2008) |arXiv:0707.0839l [hep-ph]]. 
[15] R. Easther, J. T. . Giblin and E. A. Lim, "Gravitational Waves From the End of Inflation: Computational Strategies," 

larXiv:0712.299T1 [astro-ph] . 

[16] L. R. Price and X. Siemens, "Stochastic Backgrounds of Gravitational Waves from Cosmological Sources: arXiv:0805.3570 

[astro-ph] . 

[17] S. Weinberg, Gravitation and Cosmology, John Wiley, 1972. 

[18] J. Garcia-Bellido, "Preheating the universe in hybrid inflation," arXiv:hep-ph/9804205 
[19] K. Enqvist, S. Nurmi and G. I. Rigopoulos. larXiv:0807.0382l [astro-ph]. 

[20] J. Garcia-Bellido, M. Garcia Perez and A. Gonzalez- Arroyo, "Symmetry breaking and false vacuum decay after hybrid 

inflation," Phys. Rev. D 67, 103501 (2003) [arXiv:hep-ph /0208228] . 
[21] J. Garcia-Bellido, A. D. Linde and D. Wands, "Density perturbations and black hole formation in hybrid inflation," Phys. 

Rev. D 54, 6040 (1996) [arXiv:astro-ph/9605094] . 
[22] J. Garcia-Bellido and A. D. Linde, "Preheating in hybrid inflation," Phys. Rev. D 57, 6075 (1998) [arXiv:hep-ph/9711360] . 
[23] E. J. Copeland, S. Pascoli and A. Rajantie, "Dynamics of tachyonic preheating after hybrid inflation," Phys. Rev. D 65, 

103517 (2002) [arXiv:hep-ph/0202031| . 
[24] M. Desroche, G. N. Fel der, J. M. Kratochvil and A. Linde, "Preheating in new inflation," Phys. Rev. D 71, 103516 (2005) 

|arXiv:hep-th/0501080| . 

[25] D. I. Podolsky, G. N. Felder, L. Kofman and M. Peloso, Phys. Rev. D 73, 023501 (2006) arXiv:hep-ph/0507096 . J. F. Du- 
faux, G. N. Fe lder, L. Kofman, M. P eloso and D. Podolsky, JCAP 0607, 006 (2006) arXiv:hep-ph/0602144 . 

[26] A. Buonanno, |arXiv:gr-q c/0303085. A. Buonanno, G. Sigl, G. G. Raffelt, H. T. Janka and E. Muller, Phys. Rev. D 72, 
084001 (2005) [arXiv:astro-ph/041227"7] . L. A. Boyle and A. Buonanno, Phys. Rev. D 78, 043531 (2008) |arXiv:0708.2279l 
[astro-ph]]. 

[27] S. Dimopoulos, P. W. Graham, J. M. Hogan, M. A. Kasevich and S. Rajendran , larXiv:0806.2125l [gr-qc]. 
[28] A. J. Farmer and E. S. P hinney, Mon. Not. Roy. Astron. Soc. 346, 1197 (2003) [arXiv:astro^pIi/ 0304393 . 
[29] http://www.minigrail.nl 

[30] A. Nishizawa et al, Phys. Rev. D 77, 022002 (2008) |arXiv:0710.T944l [g r-qc]]. 

[31] L. Randall, M. Soljacic and A. H. Guth, Nucl Phys. B 472, 377 (199 6) [arXiv:hep-ph/9512439| . 
[32] E. D. Stewart, Phys. Rev. D 56, 2019 (1997) [arXiv:hep-ph/9703232|. 

[33] D. H. Lyth and E. D. Stewart Phys. Rev. Lett. 75, 201 (1995) |arXiv:hep-ph/9502417| . D. H. Lyth and E. D. Stewart, 

Phys. Rev. D 53, 1784 (1996) [arXiv:hep-ph/9510204j . 

[34] R. Easther, J. T. Giblin, E. A. Lim, W. I. Park an d E. D. Stewart, JCAP 0805, 013 (2008) [arXiv:0801. 41971 [astro-ph]]. 
[35] R. Kallosh and A. Linde, JCAP 0310, 008 (2003) [arXiv:hep-th/0306058j . 

[36] S. Kachru, R. Kallosh, A. Linde, J. M. Maldacena, L. P. McAllister and S. P. Trivedi, JCAP 0310, 013 (2003) 



20 



[37] 

[38] 
[39] 



[arXiv:hep-th/0308055| . 

N. Barnaby and J. M. Clin e, Phys. Rev. D 73, 106 012 (2006) [arXiv: astro-ph /060 148 1] . N. Barnaby and J. M. Cline, Phys. 
Rev. D 75, 086004 (2007) [arXiv: astro-ph /06 1 1 750] . 

G. Felder, Comp. Phys. Comm. 179. 8 (2008) [arXiv:0712.08"T3] 

G. Felder and I. Tkachev, Comp. Phys. Comm.' 178, 12 (2008) [arXiv:hep-ph/0011159| . 



